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ABSTRACT 


This  thesis  provides  a  database  comparison  and  analysis  of  daytime 
lower  thermospheric  (E-region)  inverted  meridional  winds.  This  inversion 
technique  was  previously  developed  (Nardi,  1991)  to  retrieve  the  true  wind  at 
altitude  from  [01]  5577  A  measurements  from  the  DE-2  Satellite's  Fabry- 
Perot  Interferometer. 

The  inversion  technique  is  briefly  discussed  in  this  study.  The  raw 
DE-2  satellite  data  is  processed  to  obtain  spectrograms.  This  data  is 
processed  further  to  obtain  brightness  and  line-of- sight  winds  and  wind 
errors.  For  daytime  5577  A  emissions,  these  line-of-sight  winds  are 
averages  over  two  emission  layers.  The  purpose  of  the  inversion  technique 
is  to  "unaverage"  these  winds  through  an  onion  peel  method  on  the 
brightnesses  and  through  a  least  squares  fit  to  chapman  layers  on  the 
volume  emission  rates. 

A  data  base  of  forty-four  daytime  orbits  was  processed  through  the 
inversion  technique.  These  orbits  were  then  compared  to  two  models.  The 
Vector  Spherical  Harmonic  (VSH)  model  uses  output  from  the  National 
Center  for  Atmospheric  Research  -  Thermospheric  General  Circulation 
Model  (NCAR-TGCM)  to  produce  thermospheric  winds.  The  Horizontal 
Wind  Model  (HWM)  uses  both  empirical  data  and  a  limited  set  of  spherical 
harmonics  to  produce  thermospheric  winds.  The  output  from  both  models 
was  compared  to  the  inverted  wind  data  to  analyze  for  accuracy  and  for 
evidence  of  gravity  wave  or  tidal  activity.  Comparison's  were  also 
performed  for  specific  latitudes  and  geomagnetic  conditions. 

Vertical  profiles  were  also  produced  comparing  the  green  line  (5577 
A)  derived  meridional  winds  below  200  km  to  the  red  line  (6300  A)  winds  at 


higher  altitudes.  Again,  possible  gravity  wave  and  tidal  influences  were 
investigated. 


ACKNOWLEDGEMENTS 


Many  people  made  this  masters  thesis  possible.  First  of  all,  thanks 
goes  to  the  U.S.  Air  Force  who  has  sponsored  me  through  this  program  and 
who  now  allows  me  to  leave  active  duty.  It  is  this  ability  to  plan  and  achieve 
long-term  strategies  that  provides  us  all  with  confidence  in  the  U.S. 
Defense  Department. 

Secondly,  a  note  of  gratitude  to  my  advisor  Dr.  Timothy  Killeen  for  his 
gmdance  and  excellent  teaching  in  this  field.  Also  thanks  to  Dr.  Alan 
Burns  for  assistance  in  producing  research  products.  A  special 
appreciation  goes  to  AFIT  students  Jerry  Davis,  Tom  Smith,  Ed  Hume, 
Steve  Quigley,  and  to  Steve  Carr  for  moral  support  during  my  last 
harrowing  days  on  active  duty. 

The  staff  here  at  U  of  M  made  the  road  to  a  degree  much  easier. 
Thanks  to  Gerry  Schmitt  for  keeping  the  computers  in  line  and  on  line,  to 
Mike  Minnock  for  his  vast  computer  knowledge,  and  to  Nelly  DeMaldonado, 
who  is  by  far  the  best  secretary  and  hardest  worker  I've  ever  met. 

Lastly  thanks  goes  long  distance  to  family  and  friends  in  Atlanta.  To 
Pat  and  Jeanne  Freeman  (&  Daniel)  and  to,  of  course,  mom. 


TABLE  OF  CONTENTS 

Abstract 

Acknowledgements . i 

L  Introduction 

1.1  General . 1 

1.2  The  Lower  Thermosphere . 2 

2.  The  DE*2Fabry-Perot  Interferometer 

2.1  Introduction:  Dynamics  Explorer  Overview . 3 

2.2  The  Fabry-Perot  Interferometer . 6 

2.3  The  Interferometer . 8 

2.4  DE  -  FPI  Viewing  Geometry . 12 

3.  Greenline  •  [OB  5577  A  Emission 

3.1  Origin  of  [OB  5577  A . 14 

3.1 .1  5577  A  Nightglow . 14 

3.1.2  5577  A  Twilight/Dayglow . 15 

3.1 .3  5577  A  Aurora . 18 

3.2  Variability  of  [01]  5577  A . 19 


11 


4.  Brii^tness  Inversion 

4.1  Onion  Peel  Inversion . 21 

4.1 .1  Inversion  Method . 21 

4.1.2  Inversion  Results . 23 

4.2  Least  Squares  Fit  Inversion . 26 

4.2.1  Inversion  Method . 26 

4.2.2  Inversion  Results . 28 

4.3  Summary  and  Comparisons . 30 


5.  Wind  Inversion 

5.1  Inversion  Algorithm . 32 

5.2  Wind  Inversion  -  Weighted/Constrained . 34 

6.  Comparison  to  Model  Winds 

6.1  The  Vector  Spherical  Harmonic  (VSH)  Model . 38 

6.2  Comparisons  with  VSH  Model  Winds . 40 

6.3  Comparisons  with  HWM  Model  Winds . 43 

7.  Altitudinal  Profiles 

7.1  Red/Green  line  profile  comparisons . 46 

8.  Discussion . 53 

Appendix  A  Thesis  sub-topic:  Impact  ionization . 55 

References . 64 


CHAPTER  1 


INTRODUCTION 


1.1  General 

This  thesis  follows  a  doctoral  thesis  that  derived  an  inversion 
technique  to  recover  lower  thermospheric  meridional  winds  from  remote 
measurements  of  atomic  oxygen  5577  A  emission  (Nardi,  1991).  The 
objective  here  is  to  assemble  a  data  base  of  these  winds  and  compare  them 
to  model  output  and  to  red  line  (6300  A)  winds  from  the  same  instnunent. 

Most  of  the  inverted  wind  data  presented  here  will  be 
da)d;ime/twilight  winds.  The  inversion  technique  is  not  necessary  for 
nighttime  measurements  since  only  one  emission  altitude  exists  at  night. 

Measuring  the  lower  thermospheric  winds  has  been  a  problem  in 
thee  5577  A  range  due  to  a  double  layer  in  the  day  time.  The  satellite 
measures  the  total  emission  from  both  layers,  but  cannot  distinguish  the 
amount  contributed  to  the  total  from  each  layer.  Nardi  (1991)  has  developed 
a  method  using  a  fit  to  chapman  layers  that  allows  for  the  first  time  the 
lower  daytime  thermospheric  winds  to  be  accurately  derived  from  satellite 
measurements. 
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The  meridional  winds  are  then  compared  to  the  output  from  the 
Vector  Spherical  Harmonic  (VSH)  Model  and  the  Horizontal  Wind  Model 
(HWM).  These  models  are  presented  with  the  inverted  meridional  winds 
for  evidence  of  gravity  wave  and  tidal  effects. 

The  inverted  meridional  winds  are  also  plotted  with  higher  altitude 
red  line  (6300  A)  winds.  The  red  line  winds  require  no  inversion  process 
and  were  measured  in  an  alternating  mode  (with  the  green  line)  by  the 
satellite  late  in  the  D3mamics  Explorer  mission. 

1.2  The  Lower  Thermosphere 

The  lower  thermosphere  and  upper  mesosphere  remain  the  least 
understood  regions  of  the  atmosphere.  These  altitudes  are  too  low  for  in* 
situ  satellite  measurements  and  too  high  for  tropospheric  type  in-situ 
measurements.  Most  information  is  gathered  through  remote 
measurements  such  as  meteor  radar  (radar  tracking  of  meteor  trails  80-100 
kin),  MST  radar  (Mesosphere,  Stratosphere,  Troposphere;  echos  VHF 
signals  off  an  optically  clear  atmosphere),  and  Incoherent  Scatter  Radar 
that  sends  UHF  pulses  into  the  ionosphere. 

The  lower  thermosphere  is  affected  by  forces  that  are  not  as 
important  at  higher  altitudes.  More  radiative  energy  is  deposited  in  the 
lower  thermosphere  due  to  a  larger  optical  depth.  Tides  and  gravity  waves 
are  also  a  major  energy  source  as  they  propagate  upward  from  lower 
altiudes  and  deposit  energy  as  they  break  in  this  region. 
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CHAPTER  2 


THE  DE-2  FABRY-PEROT  INTERFEROMETER 


2.1  Introduction;  D3rnaniics  Explorer  Overview 

The  Dynamics  Explorer  (DE)  program  was  a  National 
Aeronautics  and  Space  Administration  (NASA)  sponsored  launch  of  two 
spacecraft  to  study  the  coupling  between  the  magnetosphere,  ionosphere 
and  the  atmosphere.  The  DE  spacecraft  were  launched  during  the 
summer  of  1981  and  relayed  data  until  February  1983.  DE  project 
scientists  decided  that  two  satellites  were  necessary;  one  satellite  at  a  low 
altitude  orbit  and  the  other  at  a  high  altitude  orbit.  The  satellites  were 
launched  into  coplanar  orbits  within  common  magnetic  flux  tubes.  This 
configuration  allowed  measurements  from  the  magnetosphere  and  the 
ionosphere/atmosphere  at  the  same  time.  The  DE  mission  had  four 
objectives  that  were  not  covered  in  previous  programs  (Hoffman,  et.  al., 
1981).  They  were: 

(1)  Measurements  of  the  three-component  dc  electric  fields  with 
global  coverage. 

(2)  Measurements  of  the  velocity  vector  of  the  neutral  winds  in  the 


upper  ionosphere. 
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(3)  Measurements  of  magnetic  fields,  plasmas,  particle  distribution 
functions,  neutral  gas  parameters,  characteristic  light  emissions,  and 
plasma  waves,  along  with  neutral  wind  velocities  and  electric  fields. 

(4)  Measurements  of  field  and  plasma  parameters  at  both  high  and 
low  altitudes  within  the  same  magnetic  flex  tube. 

The  instruments  on  the  DE  spacecraft  were  grouped  into  five 
categories  of  measurement  :  thermal  ions,  thermal  neutral  particles, 
suprathermal  particles,  fields,  and  optical  remote  sensing.  The  Fabry- 
Perot  Interferometer  is  a  remote  optical  sensing  instrument  and  was  a  new 
technique  for  spaceborne  measurements  when  the  DE  spacecraft  were 
launched.  Most  other  instruments  were  adapted  from  the  Atmospheric 
Explorer  (AE)  program  or  the  International  Sun  Earth  Explorer  (ISEE) 
program. 

The  primary  investigators  for  each  instrument  are  listed  in  table  2-1 . 
These  instruments  are  able  to  measure  the  following: 

-  thermal  electrons 

-  thermal  ions 

-  thermal  neutral  particles 

-  suprathermal  particles 

-  vector  magnetic  field 

-  vector  electric  field 

•  plasma  waves 

-  auroral  optical  images 

-  Very  Low  Frequency  (VLF)  transmissions  over  Antarctica 
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TABLE  2-1 


Dynamics  Explorer  Investigators 


L.H.Brace 

Goddard  Space  Flight  Center 

Langmuir  Probe 

J.L.Burch 

Southwest  Research  Institute 

HighAtitude  plasma 
instrument 

G.R.Carignai\ 

University  of  Michigan 

Neutral  mass 

spectrometer 

C.R.  Chappell 

Marshall  Space  Flight  Center 

Retarding  ion  mass 
spectrometer 

F.V.Coroniti 

University  of  California, 

Los  Angeles 

Theoretical  Program 

LA.Frank 

University  of  Iowa 

Global  auroral  imager 

W.B.Hanson 

University  of  Texas  at  Dallas 

Ion  drift  meter 

P.B.Hay8 

University  of  Michigan 

Fabry-Perot 

interferometer 

RA-Heelis 

University  of  Texas  at  Dallas 

Ion  drift  meter 

RAHelliwell 

Stanford  University 

Controlled  and 
naturally  occurring 
wave-particle 
interactions  analysis 

R.  A  Hoffman 

Goddard  Space  Flight  Center 

Low  altitude  plasma 

analyas 

N.C.Maynard 

Goddard  Space  Flight  Center 

Vector  electric  field 

instrument 

H.G.Mayr 

Goddard  Space  Flight  Center 

Theoretical  program 

AF.  Nagy 

University  of  Michigan 

Theoretical  Program 

R.G.Roble 

National  Center  for 

Atmospheric  Research 

Theoretical  Program 

E.G.SheU«y 

Lockheed  Palo  Alto  Research 

Energetic  ion  mass 

Laboratory 

spectrometer 

S.O.Shawhan 

University  of  Iowa 

Plasma  wave 
instrument 

N.W.Spencer 

Goddard  Space  Flight  Center 

Wind  and 

Temperature 

Spectrometer 

M.Sugiura 

Goddard  Space  Flight  Center 

Magnetometers 

J.D.Winningham 

Southwest  Research  Institute 

Low  altitude  plasma 

instrument 

The  apogee  and  perigee  altitudes  were  designed  to  be  high  enough  for 
auroral  viewing,  plasmasphere  and  plasmapause  observations  near  the 
equator,  and  wave  and  particle  measurements  in  the  magnetosphere. 
These  high  altitude  measurements  were  all  made  from  the  high  altitude 
spacecraft. 

From  the  low  altitude  spacecraft  measurements  needed  to  be  under 
300  km  for  neutral  particle  and  thermospheric  wind  measurements. 
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Measurements  were  also  needed  above  1000  km  for  interaction  of 
suprathermal  ions  and  the  atmosphere.  The  actual  altitude  range  of  the 
lower  altitude  satellite  was  approximately  200  km  to  1000  km. 

The  ground  data  system  for  the  DE  mission  was  designed  to  give  the 
investigators  access  to  the  information  as  qvuckly  as  possible.  DE  scientists 
used  the  same  system  that  the  Atmospheric  Explorer  mission  used.  This 
consisted  of  a  Science  Data  Processing  System  at  Goddard  Space  Flight 
Center.  A  central  computer  facility  processed  the  majority  of  the  data  at 
Goddard  with  remote  terminals  at  the  investigator's  facilities. 

2.2  The  Fabry-Perot  Interferometer 

The  Fabry-Perot  Interferometer  (FPI)  on  the  Dsmamics  Explorer  is  a 
highly  stable  single  etalon  interferometer  with  optics  able  to  view  the 
atmosphere  at  various  tangent  heights.  A  single  Doppler  profile  (spectral 
and  spatial  scan)  gives  us  the  temperature,  meridional  wind,  and  density  of 
any  one  of  several  excited  atoms  or  ions  in  the  thermosphere.  These 
Doppler  profiles  are  measured  with  a  resolution  of  10'^  A.  Table  2-2  lists  the 
species  and  spectral  regions  of  measurement. 

One  of  the  major  goals  of  the  DE  mission  was  to  increase  our 
knowledge  of  thermospheric  wind  and  temperature  fields.  How  heat  enters 
and  leaves  (is  deposited  and  dissipated)  the  thermosphere  is  not  well 
understood  due  simply  to  lack  of  data.  We  do  know  that  heat  enters  this 
region  by  solar  UV  radiation,  magnetospheric  energy  (particle  or  joule 
heating  near  the  magnetic  poles,  and  wave  or  tidal  energy  transported  from 
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TABLE  2-2 

Emissions  measured  by  FPI 


Soedes 

Soectral  Region 

Value 

0(1  D) 

6300  A 

Upper  thermosphere  wind 
and  temperature 

0+(2p) 

7320  A 

Upper  thermosphere  ion  drift 
and  temperature 

0(lS) 

5577  A 

Mesospheric  nighttime  neutral 
winds  and  lower  thermospheric 
da3rtime  winds 

N(2D) 

5200  A 

Auroral  energy  storage 

Na(2p) 

5896  A 

Mesospheric  neutral  wind  and 
temperature  daytime 

lower  altitudes.  How  much  each  of  these  contributes  and  how  their 
contributions  fluctuate  with  changing  geophysical  conditions  is  still 
unknown. 

The  Fabry-Perot  interferometer  is  shown  in  figure  2.1.  There  are 
three  sections  to  the  instrument;  the  telescope  section,  the  mirror  section, 
and  the  interferometer  section.  The  telescope  consists  of  the  main  objective 
doublet  lens  (that  is,  two  lenses  of  differing  shapes)  and  the  main  baffle. 
The  instrument  collects  light  originating  from  25  -  225  km  (tangent  point 
altitudes)  below  the  satellite.  The  FPI  is  restricted  to  looking  along  the 
satellite's  line-of-sight.  The  light  enters  the  baffle  and  stray  light  from  the 
Earth's  or  Sun's  bright  limb  or  from  cloud  tops  on  the  day-side  is 
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eliminated  by  this  baffle.  The  lens  forms  a  vertical  slice  image  of  the  sky  by 
focusing  the  collected  light  onto  the  sky  stop. 

This  vertical  slice  then  enters  the  mirror  section.  This  section 
contains  the  horizon  scan  mirror  that  is  a  concave,  front  surfaced  mirror 
that  pivots  .  This  mirror  reflects  and  refocuses  light  from  the  sky  stop  to  the 
field  stop.  The  scan  mirror  selects  different  circular  elements  from  our 
vertical  slice  of  sky.  The  pivoted  position  of  the  mirror  determines  which 
vertical  segment  of  the  sky  is  focused  onto  the  field  stop. 

The  field  stop  is  a  circular  aperture  at  the  exit  of  the  mirror  section. 
This  aperture  blocks  unwanted  portions  of  the  image  from  entering  the 
etalon  of  the  FPI  and  allows  a  small  sector  to  pass  through. 

2.3  The  Interferometer 

The  interferometer  section  of  the  FPI  consists  of  a  collimator 
telescope  lens,  a  dielectric  blocker,  a  schott  glass  blocker,  a  filter,  the  high 
resolution  etalon,  a  focusing  telescope  lens  and  the  image  plane  detector. 
Specifics  on  the  filter  wheel  and  the  scan  mirror  are  listed  in  table  2-3.  The 
emitted  light  passes  through  the  field  stop  and  is  collimated  and  filtered  for 
the  desired  wavelength.  This  light  then  passes  through  the  etalon.  The 
etalon  is  made  of  two  fused  silica  plates  separated  and  held  parallel  by  the 
three  spacing  elements.  Each  plate  is  coated  with  several  semi-reflective 
layers.  Light  passes  through  the  first  plate  and  some  light  is  reflected  back 
from  the  second  plate  while  the  rest  passes  on  through.  This  interference 
creates  a  phase  shift  in  the  light  that  creates  a  unique  "fingerprint"  or 
spectrogram  determined  by  the  intensity  and  wavelength  of 


schematic  for  the  Fabry-Perot  Interferometer  (Hays  et.  al.  1981) 
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TABLE  2-3 


rater  WheeL  1 )  Fixed  filter:  any  one  of  eight  possible  filter 

positions 

2)  Filter  cycling:  the  filter  wheel  may  be  stepped  at 
the  following  rates: 

a)  once  per  2^  *  0.5  s  (n=  0-7)  in  fast 
(default)  mode 

b)  once  per  2*^  *  0.5  s  (n=  3-7)  in  slow 
(backup)  mode 

3)  Filter  sequencing:  the  filter  wheel  may  be 
commanded  to  alternate  between  any  two  filter 
positions 

4)  Impulsive  go  step  (automatic  filter  sequence)  if 
logic  is  lost. 

Scan  Mirror;  1)  Fixed  mirror:  any  one  of  1 6  sky  positions  and  2 

calibration  positions 

2)  Scan  mirror  cycling:  the  mirror  may  be 
stepped  once  per  (2”  *  0.25)  +  0.25  s  (n  s  0-7) 

3)  Impulsive  go  step  (automatic  scan  sequence)  if 
logic  is  lost 


the  light.  The  etalon  transforms  the  spectral  distribution  of  an  emission 
line  into  a  radial  distribution  in  the  form  of  a  circular  fringe  pattern.  This 
circular  pattern  is  what  the  image  plane  detector  measures. 

The  image  plane  detector  (IPD)  converts  the  light  intensity 
distribution  of  the  interference  fringe  into  discreet  electron  pulses.  These 
pulses  are  measured  a  segmented  anode  made  up  of  12  concentric  rings  of 
equal  area.  The  measured  output  of  each  anode  ring  corresponds  to  one  of 
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figure  2.2  The  multi-channel  detector's  measurement  of  the  fringe  pattern  produced  by  the  etalon. 
One  radial  'cut'  of  one  order  of  the  fringe  will  produce  a  spectrogram  that  gives  wind  and  temperature 
information.  The  multi-channel  detector  samples  the  entire  fringe  simultaneously  for  greater 
resolution. 
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12  channels  on  a  spectrogram.  Figure  2.2  illustrates  this  concept.  From 
these  spectrograms  we  can  extract  wind,  temperature,  and  brightness 
information. 

2.4  DE  -  FPI  Viewing  Geometry 

The  sixteen  sky  positions  of  the  horizon  scan  mirror  allow  a  scan 
from  5°  -  15“  below  the  spacecraft's  local  horizontal  (tangent  to  the  Earth). 
Figure  2.3  is  a  schematic  showing  the  viewing  geometry  of  the  FPI  and  DE 
satellite.  The  tangent  point  is  the  point  along  the  line  of  sight  view  of  the 
FPI  where  the  line  of  sight  is  tangent  to  a  shell  concentric  with  the  surface 
of  the  earth.  This  height  happens  to  be  the  altitude  where  the  contribution 
to  any  FPI  measurement  is  at  a  maximum.  This  is  because  the  path  length 
of  the  line-of-sight  is  longest  through  a  layer  centered  at  the  tangent  point. 

The  tangent  heights  range  from  ~25  -  225  kilometers  below  the 
satellite  altitude.  The  instrument  axis  is  an  imaginary  line  passing 
through  all  of  the  optical  surfaces  of  the  instrument.  The  instrument  axis 
lies  10.294“  below  the  local  horizontal  for  the  majority  of  the  DE-2  mission. 

Finally,  the  view  axis  is  the  same  as  the  line-of-sight  and  varies  from 
5“  -  15“  below  the  local  horizontal.  This  is  the  imaginary  line  that  passes 
through  the  centroid  field-of-view  for  any  viewing  position. 


StteOiie  Measuronent  Geomeoy:  OE2-FPI 


figura  2.3  OE2-FPI  viewing  geometry.  This  schematic  shows  the  instrument  axis,  the  range  of  the 
horizon  scan  and  the  varying  field-of-view  for  the  OE2>FPI. 
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CHAPTERS 


GREENLINE  -  OI 5577  A  EMISSION 

3.1  Origin  of[Ori  5577  A 

The  5577  A  emission  originates  from  the  ->^0  electric  quadrapole 
transition  of  atomic  oxygen: 

0(lS)  ->  0(lD)  +  hv  (5577  A)  (3.1 ) 

3.1.1  5577  A  Nightglow 

The  01  5577A  nightglow  has  two  layers;  one  at  97  ±  2  km  and  an 
upper  layer  at  approximately  240  km.  The  upper  layer  is  one  order  of 
magnitude  dimmer  than  the  lower  layer  (Cogger  and  Anger,  1973).  The 
peak  volume  emission  rate  of  the  upper  layer  is  3  -  4  photons  sec'^  cm-3 
compared  to  50  -  360  photons  sec'^  cm*3  for  the  lower  layer. 

For  the  lower  layer  (<  120  km)  the  product  of  0(^8)  is  thought  now  to 
be  due  to  the  Barth  mechanism; 


0  +  0  +  M  O2*  +  M 

(3.2) 

02*+  0  ->02  +  0(1S) 

(3.3) 

where  the  O2*  is  most  likely  02(c^Xu')  (Bates,  1988). 
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For  the  upper  layer  dissociative  recombination  is  the  source  of  0(lS)  green 
line  airglow: 

O2+  +  e  ->  O2*  0(1S)  +  0  ( 3.4 ) 

The  O2*  can  be  produced  in  one  of  two  ways: 

O2+  +  e  02(cllu-)  ->  0(1S)  +  0(1D)  +  0.83  ev  ( 3.5 ) 

O2+  +  e  ^  02(3ng)  ->  0(1S)  +  0(3P)  +  2.79  ev  ( 3.6 ) 

The  question  as  to  which  of  these  branches  is  primarily  responsible 
for  the  production  of  0(1  S)  is  still  under  study.  Branching  ratio  studies 
were  performed  (Hernandez,  1971)  using  ground  based  Fabry-Perot 
interferometer  data.  These  studies  suggested  equation  (  3.6  ).  The  problem 
with  this  study,  however,  is  that  the  upper  layer  of  5577A  emission  is 
masked  by  the  brighter  lower  layer  emission  at  97  km. 

Killeen  and  Hays  ( 1981  )  used  DE2-FPI  data  from  the  DE2  satellite  to 
obtain  accurate  branching  ratio  values.  Their  results  favored  equation 
(  3.5  ).  The  recombination  reaction  that  leads  to  the  production  of  0(1  S)  is 
dependent  on  the  electron  temperature  (  Zipf,  1988  )  and  density  (Abreu  et 
al.,  1983)  as  well  as  the  vibrational  distribution  of  02"*^  (Guberman,  1979; 
Zipf,  1979,  1980).  Inconsistences  arise  in  the  populations  and  collisional 
frequencies  of  02'*'  • 

3.1.2  5577A  Twilight^Dayglow 

Sources  for  the  5577A  dayglow  are  the  same  as  for  the  nighttime 
sources  plus  others  as  summarized  below: 
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Nlghtgtow; 

0  +  0  +  M  -4  O2*  +  M 

(3.2) 

0 

+ 

0 

->02  +  0(1S) 

(3.3) 

and 

02'*’  +  e  ->  O2*  0(1S)  +  0 

(3.4) 

Others: 

Direct  Photodissociation: 

02  +  hv  ->  0(1S)  +  0 

(3.7) 

Photoelectron  impact  excitation: 

0  +  e*  -4  0(1S)  +  e* 

(3.8) 

Quenching  of  excited  N2: 

N2(A3Iu+)  +  0(3P) 

0(1S)  +  N2(Xl2:g+) 

(3.9) 

Atom-ion  interchange: 

N  +  O2+  ->  0(1S)  +  NO+ 

(3.10) 

All  of  these  sources  Oie  significant  in  twilight  airglow.  Which 
source  is  dominant  depends  upon  the  altitude. 

The  lower  layer  at  ~100  km  is  still  present  in  the  twilight  and  daytime 
airglow,  but  is  somewhat  thicker  than  at  night.  This  is  due  to  reaction 
(3.7),  the  production  of  0(1  S)  through  direct  dissociation  of  molecular 
oxygen.  This  is  the  dominant  reaction  above  260  km  and  significant  above 
190  km. 

There  is  also  an  upper  layer  in  twilight,  but  the  peak  is  lower  than  at 
night.  In  fact,  the  upper  layer  is  highest  at  night  (~240  km)  and  descends 
in  altitude  with  decreasing  solar  zenith  angle  as  da3rtime  conditions 
develop.  The  twilight  altitude  of  the  upper  layer  is  ~200*220  km  and  the 
daytime  peak  is  at  -160-170  km.  The  intensity  of  the  upper  layer,  however, 
increases  from  3-4  photons  sec’l  cm-3  at  night  to  -10-60  photons  sec'l  cm-3 
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figure  3.1  Volume  emission  rates  from  orbits  7461  and  6165,  Two  layers  of  emission  can  bo  seen, 
one  at  97  km  and  one  at  -1 60  Km.  The  upper  layer  decreases  in  intensity  and  raises  in  altitude  with 
increasing  solar  zenith  angle. 
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for  twilight  to  ~1 20-200  photons  sec^  cm*3  in  dayglow.  Figure  3.1  shows  two 
surface  plots  of  volume  emission  rates  as  a  function  of  angle  along  track  (or 
latitude)  and  altitude.  From  these  plots  it  is  evident  that  there  are  two 
layers  of  emission  and  that  the  upper  layer  increases  in  intensity  as 
daytime  conditions  develop  and  decreases  with  increasing  solar  zenith 
angle. 


These  plots  do  not  show  the  change  in  altitude  of  the  upper  layer  very 
well  because  they  have  constrained  the  volume  emission  rates  to  two 
specific  altitudes,  97  and  160  km.  These  constrained  volume  emission  rates 
will  be  discussed  in  detail  in  the  next  chapter. 

3.1.3  5577A  Aurora 

The  emission  of  [01]  5577A  during  aurora  is  extremely  variable  in 
intensity  and  altitude.  Because  of  this,  its  sources  are  difficult  to 
determine.  Whereas  the  nightglow  intensity  is  ~60-250  Rayleighs  (lower 
layer)  emd  the  dayglow  intensity  is  ~200-1000  Rayleighs,  the  auroral  5577 A 
emission  is  -500-7000  Rayleighs  with  a  peak  at  105  km  (Sharp  et  al.,  1979; 
Solomon  1987). 

Much  study  has  been  accomplished  to  determine  which  reaction  is 
the  dominant  source  of  the  emission  with  significant  uncertainty 
remaining.  Reaction  (3.9)  is  now  accepted  as  the  major  producer  of  OC^S) 
in  the  aurora  (Sharp  et  al.,1979;  McDade  and  LLewellyn,  1984b,1985).  Other 
reactions  have  been  considered,  but  are  thought  to  be  only  minor  sources  of 
0(1S). 
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One  other  possible  reaction  still  iinder  study  is: 

O2*  +  0(3P)  ^02  +  0(1S)  ( 3.11 ) 

This  source  was  suggested  by  Yau  and  Shepherd  (1 979)  and  could  be 
a  major  source  below  105  km.  Here,  the  O2*  is  probably  02(cllu‘)  from 
auroral  high  energy  electrons  (Solheim  and  LLewell5m,  1979). 

3.2  Variability  of[OI]  5577  A 

The  latitudinal,  seasonal,  diurnal  and  solar-cycle  variations  in  the 
[01]  5577  A  are  quite  pronounced.  Before  any  analysis  of  lower 
thermospheric  winds  can  be  performed,  a  basic  understanding  of  these 
variations  is  necessary. 

The  seasonal  variation  appears  to  be  semiannual  in  nature.  The 
maximum  emissions  were  found  to  be  just  after  equinoxes  (April  •  May, 
October  -  November)  and  minimum  emissions  were  just  after  the  solstices 
(January  -  February,  July  -  August).  Smith  and  Steiger  (1968)  foimd  a 
variation  factor  of  6  in  from  ~20°  to  32®  N  latitude. 

Subsequent  ground  based  and  satellite  studies  have  confirmed  that 
this  semiannual  variation  is  somewhat  dependent  on  latitude  with 
minimum  variations  at  the  geographic  equator  and  maximvim  variations 
from  40®  to  60®  in  both  hemispheres. 

The  diurnal  variation  is  semidiurnal  just  as  the  seasonal  variation 
was  semiannual.  Studies  have  shown  a  maximum  increase  of  50-100%  in 
the  5577  A  intensity  in  this  semidiurnal  variation  caused  by  the 
semidiurnal  solar  tide.  The  maximum  ranges  in  local  solar  time 
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an)rwhere  from  22  -  06  LST.  The  exact  time  depends  on  the  month  and 
latitude  but  still  shifts  from  year  to  year.  It  is  because  of  this  shifting  that 
long  term  studies  did  not  yield  any  clear  patterns  since  the  maximums 
were  averaged  out  and  produced  no  clear  daily  oscillation. 

The  solar  cycle ! geomagnetic  variation  is  most  apparent  from  30°  to 
60°  latitude.  This  variation  does  seem  to  be  dependent  on  solar  cycle/solar 
sunspot  number.  The  variation  is  also  dependent  on  the  orientation  of  the 
interplanetary  magnetic  field  (IMF)  with  an  increase  of  nightglow 
measurements  under  positive  IMF  conditions  and  decreased  intensities 
under  negative  IMF  conditions. 

Small  Scale  Variations  are  evident  in  the  measurements  as  a  result 
of  gravity  waves  and  tides.  Waves  are  an3rwhere  from  6-250  km  with 
transverse  lengths  as  high  as  hundreds  of  kilometers. 
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CHAPTER4 


BRIGHTNESS  INVERSION 


4.1  Onion  Peel  Inversion 

4.1 .1  Inversion  Method 

The  idea  behind  the  'onion  peel'  method  is  to  'peel'  away  brightness 
measurements  from  upper  atmospheric  layers  successively  to  reveal 
contributions  from  lower  layers.  We  assume  or  know  nothing  about  the 
amount  of  contribution  these  lower  layers  might  yield  before  performing  the 
inversion.  Altitude-scanned  limb  brightness  measurements  are  used  from 
the  DE2-FPI.  An  inversion  algorithm  solves  for  the  inverted  or  'onion 
peeled'  volume  emission  rate  by  relating  the  height  integrated  volume 
emission  rate  to  the  limb-scanned  brightness  measurements. 

In  terms  of  an  equation: 

B(zt)  =  I  K(2,zt)  Ti(z)  dz  ( 4.1 ) 

z 


B(zt)  =  Surface  brightness  at  tangent  height  zt 

T]  =  Volume  emission  rate  profile  (function  of  altitude  ,  z) 

K(z,zt)  =  Brightness  Inversion  kernel: 
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K(zt,z)  =  X  Ka/fi  (zt,z)  ( 4.2  ) 

M 

Ka/FI  =  Altitude  broadening  kernels  for  the  individual  ring  anodes  , 
channels  1-12  from  the  FPI  image  plane  detector.  These  kernels  are 
functions  of  the  field-of-view  weighting  functions  from  raytrace 
simulations.  These  raytrace  simulations  were  necessary  to  correct  for 
distortions  produced  by  the  FPI’s  optics.  The  raytrace  takes  an  individual 
ray  of  light  from  its  emission  source  in  the  atmosphere  and  follows  it 
through  the  FPI  optical  components  to  where  it  strikes  the  optical  sensor  of 
the  FPI  (see  figure  2.1 ). 

To  actually  perform  calculations  using  equation  4.1,  it  must  first  be 
represented  in  matrix  form; 

[B]  s  [K][n]  ( 4.3 ) 

Here,  [B]  and  [r|]  are  column  vectors  (Nxl)  and  [K]  is  an  NxN  matrix. 
The  next  step  is  to  invert  the  [K]  matrix  to  get  the  solution  [t|]: 

[Til  =  [K]-i  [B]  ( 4.4 ) 

This,  however,  is  best  done  by  using  sing^ar  value  decomposition 
(SVD).  SVD  rewrites  or  'decomposes'  the  [K]  matrix  into  three  other 
matrices: 

[K]  =  [U][W][V]T  (4.5) 

Where  [U]  is  the  eigenvector  matrix  of  [K]  and  [V]  is  the  eigenvector 
T 

matrix  of  [K]  and  [U]  and  [V]  are  orthogonal.  Performing  the  inversion 
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this  way  helps  to  deal  with  singularities  or  near-singularities  in  the  kernel 
matrix  [K]  by  isolating  the  problems  in  the  eigenvalue  matrix  [W].  Since 
[U],[V],  and  [W]  are  all  either  orthogonal  or  diagonal,  they  are  much  easier 
to  manipulate  than  the  [W]  matrix.  The  final  solution  can  now  be  written 
as: 

(Til  =  [K]‘  [B]  =  {(U][W-][V]T}-'  [B] 

=  [V]  •  [diag(l/[Wj ))]  ■  (rW^  •  [B]  )  ( 4.6 ) 


4,1.2  Inversion  Results 

Figure  4,1  shows  the  unconstrained  brightness  inversion  profile 
(plusses).  The  fit  is  close  to  the  model  brightness  everywhere  except  at 
large  vertical  gradients  (the  E-region  layer).  Note  that  in  the  upper  daytime 
plot  the  F-region  layer  is  resolved  well  with  this  ‘onion  peel’  inversion.  The 
poor  resolution  of  the  E-region  layer  is  also  evident  in  the  bottom  panel 
where  there  is  no  upper  F-region  emission  layer. 

Figure  4.2  shows  a  surface  plot  of  the  unconstrained  inversion  at  the 
top.  The  bottom  surface  plot  is  a  2-layer  chapman  fit  (to  be  discussed  next 
section).  The  E-region  once  again  is  not  as  well  resolved  in  the  upper  plot 
as  it  is  in  the  lower  chapman  layer  plot.  The  reason  for  this  is  that  the 
chapman  layer  inversion  constrains  the  emission  layers  to  specific 
altitudes  with  specific  half-thicknesses  using  a  least  squares  algorithm. 
This  is  the  topic  of  the  next  section. 
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Volume  Emission  Rate  [ph/cc*s] 


Volume  Emission  Rate  [ph/cc-s] 


figure  4.1 
'Onion  Peel' 
brightness  inversion 
for  daytime  and 
nighttime  conditions. 
The  line  is  a  plot  of  a 
model  volume 
emission  rate,  the 
boxes  a  simulated 
limb  brightness 
profile,  and  plusses 
are  the  inverted 
volume  emission 
rate.  The  inverted 
inversion  rate  fits  the 
model  well  where  no 
large  vertical 
gradients  exist.  The 
E  region  peak, 
however,  is  not 
handled  as  well  by  the 
inverted  volume 
emission  rate  where 
there  is  a  large 
vertical  gradient  due 
to  near  singularities 
caused  by  the  large 
field  of  view  of  the 
DE2-FPI. _ 

□  simulated  limb  teightness 
invened  volume  emission  rate 

—  model  volume  emission  rate 
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2000 


figur*  4.2  Brightnass  inversions  for  DE2  orbit  6832.  The  top  surface  plot  shows  the  unconstrained 
‘onion-peer  inversion  whiie  the  bottom  plot  shows  the  2-layer  chapman  fit  inversion. 
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4.2  Least  Squares  Fit  Inversion 


4,2.1  Inversion  Method 


An  alternative  to  performing  an  \mconstrained  inversion  is  to  pre¬ 
determine  the  shape  of  the  volume  emission  rate  profile.  The  layers  of 
emission  are  known  as  Chapman  layers  with  assumed  peak  heights  and 
half- thicknesses. 

From  a  known  Chapman  emission  layer  we  produce  simulated 
limb-scan  brightness  profiles.  These  simulated  profiles  are  then  fitted  to 
DE2-FPI  limb-scanned  brightness  measurements  to  determine  the  best 
Chapman  layer  ciiaracteristics.  The  volume  emission  rate  profile  is 
obtained  bv  performing  a  scaled  summation  of  the  Chapman  emission  rate 
layers. 

This  chapter  will  address  a  two  layer  Chapman  layer,  but  as  many 
layers  as  desired  can  be  used  (i.e.  four, six,...).  The  two  Chapman  layer  has 
the  following  characteristics: 


-Region _ 

E-Region  (night) 
E-Region  (day) 
F-Region  (night) 
F-Region  (day) 


Peak  Emission 

97  km 
97  km 
200+  km 
160  km 


Half-thickness 

10  km 
-26  km 
Variable 
Variable 


TABLE  4-1 
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The  mathematics  behind  the  inversion  begins  with  the  formula  for 
the  model  Chapman  function  emission  layers: 


Ti^..(z)*Y.  exp 


H,,  »»,  H,, 


1  =  1.2 


(4.7) 


Tici(z)  =  Volume  emission  rate  for  Chapman  layer  i 
Hsi  =  Fixed  bottom-side  scale  height 
Hxi  =  Fixed  top-side  scale  height 
zi  s  Fixed  height  of  peaks  emission 
Yi  =  Arbitrary  magnitude 


The  simulated  limb  brightness  profile,  Bci(zt)*  1=1.2,  that  matches 
this  model  Chapman  function  is: 


Bci(zt)  =  /  K(z,zt)  Tlci(z)  dz,  i  =  1,2  ( 4.8  ) 

FOV 

Here  K(z,zt)  is  the  same  brightness  kernel  as  discussed  in  the  last 
section.  The  two  simulated  brightness  profiles  Bci(zt)  and  Bc2(zt)  are  fitted 
using  constants  Xi  and  X2  to  match  as  closely  as  possible  (least  squares 
algorithm)  the  measured  DE2-FPI  brightness  profile  B(zt)  as  in: 


[Be]  •  [X]  =  [B] 


(4.9) 
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The  least  squares  method  is  used  to  obtain  the  best  solution  for  the 
modulating  constants,  Xi,  i=l,2.  Once  these  constants  are  known,  the 
inverted  volume  emission  rate  profile  is  found  by  calculating  the  weighted 
sum  of  the  original  Chapman  layers  from: 

Tl(z)  =  Xi  Tlcl(z)  +  X2  T1c2(z)  ( 4.10  ) 

In  other  words,  if  the  Xi  layer  had  a  contribution  of  1  and  the  X2 
layer  had  a  contribution  of  2,  then  the  t1c2(2)  volume  emission  rate  layer 
would  be  twice  as  large  as  the  TidCz)  emission  rate  layer.  Since  the  X’s  are 
directly  dependent  on  the  Chapman  layer  form,  it  is  evident  that  the 
Chapman  layer  chosen  has  a  direct  effect  on  the  shape  of  the  resulting 
volume  emission  rate. 


4.2.2  Inversion  Results 

Figure  4.3  shows  the  two  brightness  inversions  (imconstrained  and 
2-layer  chapman  fit)  for  orbit  6565.  Here  it  is  easier  to  see  the  poorer 
resolution  of  the  E-region  by  the  unconstrained  inversion  in  the  top  surface 
plot.  The  F  and  E  regions  are  distinguishable,  but  not  to  the  extent  of  the 
chapman  inversion  shown  on  the  bottom  plot.  Once  again  with  both  plots 
the  F  region  emission  decreases  in  magnitude  as  twilight  approaches.  In 
the  top  plot,  however,  the  lift  in  altitude  can  also  be  seen.  The  chapman 
inversion  cannot  show  this  because  we  have  predetermined  the  two 
emission  layer’s  altitudes  and  these  altitudes  cannot  change  without 
running  another  inversion  with  new  emission  layer  altitudes. 
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figure  4.3  Brightness  inversions  for  DE2  orbit  6565.  The  unconstrained  inversion  is  at  the  top  while 
the  2-(ayer  Chapman  fit  is  at  the  bottom.  The  better  resolution  of  the  E-region  is  evident  in  the 
Chapman  layer  fit. 


4.3  Summary  and  Comparisons 


The  following  table  summarizes  the  two  inversion  techniques: 

INVERSION  TECHNIQUES 


rMJCiimitlJCi 


2-LAYER  4-LAYER 


Daytime  F-Region 


Da3rtime  E-Region 


TABLE  4-2 

The  night-time  inversion  is  not  resdly  needed  since  the  only  emission 
layer  is  the  E-region  layer.  No  upper  F-region  layer  exist  to  add  to  the  line- 
of-sight  winds,  so  the  actual  line-of-sight  winds  are  the  actual  E-region 
winds. 

The  advantages  of  the  chapman  layer  inversion  is  that  once  again  the 
E-region  is  better  resolved.  But  there  are  several  disadvantages.  The  main 
problem  is  that  it  is  difficidt  to  know  the  exact  height  and  half-thicknesses 
of  the  emission  layers.  Another  problem  occurs  when  the  satellite  altitude 
varies  more  than  +/-  30  km.  The  assumed  chapman  emission  layers  are 
calculated  with  a  fixed  satellite  altitude  in  mind.  When  the  actual  altitude 
varies  from  this  altitude  too  much  the  accuracy  of  the  fit  is  compromised. 
Of  course  we  could  just  change  the  reference  chapman  layer  brightness 
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profiles  with  changing  satellite  altitude,  but  this  would  introduce  much 
more  computational  time. 

The  unconstrained  inversion,  on  the  other  hand,  is  less 
computationally  extensive  and  more  desirable  even  though  the  E-region  is 
not  resolved  as  well.  This  inversion  also  calculates  the  brightnesses  at  the 
tangent  heights  which  is  exactly  what  is  required  for  the  wind  kernel 
evaluations.  Lastly,  the  unconstrained  inversion  gives  maximum  freedom 
to  the  shape  of  the  emission  profile.  That  is,  it  allows  the  altitude  shift  of 
the  emission  layers  from  day  to  night. 

Because  of  these  three  advantages,  the  unconstrained  ‘onion  peel’ 
inversion  is  used  to  calculate  the  wind  inversion  kernels. 
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CHAPTERS 


WIND  INVERSION 


5.1  Inversion  Algorithm 

The  wind  inversion  algorithm  is  a  weighted/constrained  inversion. 
The  volume  emission  rate  profiles  used  to  calculate  the  wind  inversion 
algorithm  are  acquired  from  a  previous  onion  peel  inversion  of  brightness 
measurements  as  discussed  in  the  previous  chapter. 

The  line-of-sight  (LOS)  winds  are  averaged  winds  over  the  viewing 
path  of  the  DE  satellite.  Figure  5.1  illustrates  this  concept.  To  get  the 
actual  wind  at  each  altitude,  another  inversion  process  must  be  applied. 

The  first  step  is  to  write  an  equation  that  relates  the  line-of-sight 
wind  to  the  raw  data  of  the  DE2-FPI.  This  has  been  done  through  a 
simplified  variation  of  a  four  parameter  least  squares  fit  algorithm.  After 
this,  a  relationship  between  the  line-of-sight  winds  and  the  actual  winds  at 
height  z  must  be  determined.  This  involves  relating  changes  in  the  line-of- 
sight  wind  to  changes  in  the  response  of  the  channels  of  the  FPI  through  a 
least  squares  method.  The  changes  in  the  response  of  the  channels  of  the 
FPI  include  the  volume  emission  rate  profiles  previously  obtained  (and 


LINE-OF-SIGHT  WIND  MEASUREMENT:  DE-FPI 
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produced  each  unaveraged  doppler  shift. 
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associated  brightness  kernels).  The  final  line-of-sight  wind  equation  can  be 
written: 


vios  (zt)  =  J  K'(z,zt)  V  (z)  dz  +  viosO  (zt)  (  5.1 ) 

where: 

vios  (zt)  =  Line-of-sight  wind  at  tangent  height  z 

K'(z,zt)  =  Wind  inversion  kernel  (includes  altitude  averaging 
effect) 

viosO  (zt)  ~  Apparent  winds  -  these  are  artificial  winds  produced 
from  gradients  in  the  volume  emission  rate  profile 

The  wind  kernel  K'  also  includes  a  wind  offset  factor  of  330  m/s  for 
5577A  emission.  This  offset  is  due  to  technical  factors  related  to  how  the 
FPI  measures  the  brightness  data. 


5.2  Wind  Inversion  -  weighted/constrained 


This  inversion  is  weighted  by  a  covariance  matrix  Syios  as  shown: 


v  =  (yH  +  KTSvlos-^K)-!  (KTSvios-lvios) 


(5.2) 
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where: 

V  s  actual  wind 

Y  =  A  weighting  factor  whose  optimal  value  is  0.003  .  This 

value  was  determined  by  testing  various  combinations  of  y 
and  the  Nth  order  difference  matrix  H. 

H  s  An  nth  order  difference  matrix.  This  difference  matrix 
limits  the  munber  of  inflection  points  in  the  inverted  wind 
profile  and  here  n=l.  The  higher  the  value  of  n,  the  less 
constraint  on  the  solution. 

K  s=  The  wind  inversion  kernel 

Svios  =  The  covariance  matrix  consisting  of  the  line-of-sight 
wind  errors 

vios  =  The  line-of-sight  winds 

The  accuracy  of  the  inverted  winds  using  this  weighted  constrained 
procedure  are  +/- 16  m/s  with  a  spatial  resolution  of  +/- 10  km.  Figure  5.2 
shows  inverted  winds  for  three  profiles  of  DE2  orbit.  The  fine  wave 
structure  is  evident  in  the  inverted  wind,  particularly  at  lower  altitudes  (~ 
100 -150  km). 

Figure  5.3  shows  surface  plots  of  the  inverted  wind  for  orbits  6565, 
6725,  and  6745.  A  definite  fine  wave  structure  can  be  seen  both  altitudinally 
and  latitudinally.  These  are  probably  the  results  of  gravity  waves  and/or 
tides  propagating  from  lower  altitudes. 
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Wei^ted/C<MistraiDed  Wind  Inversioii 
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figure  5.2  Three  consecutive  profiles  of  DE2  orbit  6565  showing  the  line-of-sight  winds  and  the 
inverted  winds  on  the  right  and  associated  wind  kernels  on  the  left. 
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M^dional  Winds:  Weighted  Constrained  Inversion 


figure  5.3  Three  surface  plots  of  the  inverted  meridional  winds  from  orbits  6565,  6725,  and  6745. 
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CHAPTER6 


COMPARISONS  TO  MODEL  WINDS 


6.1  The  Vector  Spherical  Harmonic  (VSH)  Model 

The  Vector  Spherical  Harmonic  (VSH)  model  is  a  semi-empirical 
model  that  uses  output  from  the  NCAR-TIGCM.  The  NCAR-TIGCM 
(National  Center  for  Atmospheric  Research's  Thermospheric  Ionospheric 
General  Circulation  Model)  solves  the  three  dimensional  time  dependent 
equations  that  describe  the  thermosphere. 

The  VSH  model  uses  TIGCM  derived  coefficient  files.  These 
coefficients  are  actually  weighting  factors  created  from  operations  on  the 
TIGCM  wind,  temperature,  and  composition  outputs.  This  operation  is  a 
vector  spherical  expansion  on  vrind  fields  and  a  scalar  spherical  expansion 
on  the  temperature  and  composition. 

These  coefficient  files  reside  in  a  coefficient  library  and  include  a 
variety  of  conditions  (seasonal  variation,  magnetic  activity,  and  solar 
activity).  These  files  total  to  18  nms  of  the  TIGCM  ready  for  use  by  the  VSH 

model. 
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INVERTED  MERIDIONAL  WIND  VS  VSH  WIND  •  110  km 
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figure  6.1  TTiree  Oi-bits  of  DE2  showing  the  surface  plot  to  the  left  and  an  altitude  slice  at  110  km 
comparing  the  inverted  meridional  wind  (boxes)  to  the  VSH  model  wind  (plusses). 
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6.2  Comparisons  with  VSH  Model  Winds 

Figure  6.1  shows  altitude  slices  at  110  km  for  three  DE2  orbits 
comparing  the  inverted  meridional  winds  to  the  VSH  winds.  It  is 
immediately  evident  that  there  is  a  fine  wave  structure  to  the  meridional 
winds  that  the  VSH  model  is  not  capable  of  producing.  These  waves  are  the 
probably  the  result  of  gravity  waves  or  tidal  motions.  The  NCAR-TIGCM 
has  been  modified  to  take  into  account  the  effects  of  upward  propagating 
tides  in  the  lower  thermosphere  (  Fesen  et  al.,  1986).  The  effect  of  this 
modification,  however,  is  limited  in  the  VSH  model  and  shows  up  as  a 
semidiurnal  effect  on  the  meridional  wind  output.  Gravity  waves  are  not 
well  understood  and  still  very  difficult  to  predict  or  model. 

Over  forty  DE2  orbits  were  processed  to  provide  inverted  meridional 
green  line  winds.  These  winds  were  binned  and  averaged  by  latitude  and 
ap  indices  and  are  shown  in  figure  6.2.  The  low  ap  averages  show  the  best 
correlation  to  the  VSH  model  meridional  winds,  especially  between  -45®  and 
45®  latitude.  At  higher  latitudes,  auroral  effects  cause  much  greater 
perturbations  in  the  inverted  winds. 

The  medium  ap  averages  also  show  a  good  correlation  with  the 
model  output,  except  at  high  latitudes.  The  high  ap  VSH  output,however, 
shows  a  fairly  poor  match  to  the  inverted  winds  at  all  latitudes. 

These  results  are  what  might  be  expected  given  previous  studies  of 
the  dynamics  of  the  lower  thermosphere.  Cogger  et  al.  (1985)  studied  the 
5577A  emission  over  Calgary,  Canada  (51°  N,  114°  W,  geographic)  using  75 
nights  of  data.  He  found  a  highly  variable  wind  field  with  maximum 
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figure  6.2  Comparisons  of  VSH 
model  meridional  wind  (plusses) 
to  the  averaged  inverted 
meridional  wind  (boxes)  at  110 
Km.  The  Low  ap  at  mid  and  low 
latitudes  shows  the  best 
correlation  while  the  high  ap  at 
high  latitudes  shows  the  worst. 
The  averages  are  binned  over 
10®  latitude  and  represent 
approximately  forty  DE2 
orbits. 


□  •  Averaged  inverted 

Meridional  Wind 


-|-  •  VSH  Meridional  Wind 
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figure  6.3  Averaged  low  ap 
meridional  winds  vs.  VSH 
meridional  winds  at  110  km. 

The  top  plot  shows  an  average 
for  ten  orbits  mostly  at  mid  and 
low  latitudes.  The  bottom  plot 
includes  twelve  more  orbits  at 
high  latitudes  under  auroral 
conditions.  In  the  bottom  plot, 
auroral  conditions  cause  much 
more  variation  in  the 
meridional  wind  that  the  VSH 
model  cannot  match. 
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speeds  near  200  ms'i  during  disturbed  geomagnetic  periods  At  times, 
stable  wind  patterns  were  observed,  at  other  times,  little  organization  was 
present. 

Johnson  and  Luhmann  (1985)  studied  MST  (Mesospheric, 
Stratospheric,  Tropospheric)  radar  data  from  Poker  Flat,  Alaska  (65°N). 
They  found  the  thermospheric  winds  to  be  poorly  ordered  and  highly 
variable  in  relation  to  geomagnetic  activity.  The  results  presented  here 
seem  to  be  in  agreement  with  these  findings. 

Figure  6.3  shows  two  versions  of  averaged  low  ap  inverted 
meridional  wind.  The  top  plot  shows  a  mixture  of  all  conditions  while  the 
bottom  plot  includes  auroral  conditions.  The  variability  is  greater  for 
auroral  conditions  as  expected,  although  again  the  best  correlation  to  the 
model  output  occurs  at  low  and  mid  latitudes. 

6.3  Comparison  to  the  HWM  Model 

The  inverted  meridional  winds  are  also  compared  to  the  Horizontal 
Wind  Model  (HWM90).  The  HWM90  is  a  revised  version  of  the  HWM87 
empirical  model  (Hedin  et  al.,  1988)  that  provides  upper  thermospheric 
winds  based  on  DE2  and  AE-E  satellite  data.  The  newer  version 
incorporates  ground-based  incoherent  scatter  radar  and  Fabry-Perot  data  to 
extend  the  lower  range  of  the  model  from  220  km  to  100  km.  A  limited  set  of 
spherical  harmonic  expansions  are  used  in  this  model.  The  HWM  provides 
winds  in  a  manner  similar  to  the  MSIS-86  model's  output  of  temperature, 
density,  and  composition. 

The  model  agrees  well  with  various  data  sets,  especially  at  low  and 
mid  latitudes  (Hedin  et  al.,  1991).  Solar  activity  variations  are  also  included 
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in  the  model,  although  they  were  found  to  be  small.  The  model  describes 
diurnal  variations  in  the  upper  thermosphere  and  semidiurnal  variations 
in  the  lower  thermosphere. 

Figure  6.4  shows  three  orbits  of  DE2-FPI  inverted  green  line  data. 
The  surface  plots  are  shown  on  the  left  with  an  altitude  slice  plotted  at  the 
right.  Also  plotted  are  the  VSH  model  meridional  winds  (plusses)  and  the 
HWM  model  meridional  winds  (triangles).  Consistent  with  previous 
findings  of  the  HWM  model,  the  best  fit  to  the  inverted  winds  occur  at  mid 
and  low  latitudes.  The  VSH  and  HWM  models  do  not  disagree  with  each 
other  to  any  great  extent,  but  both  clearly  miss  the  effects  of  gravity  waves 
and  tides  sen  in  the  inverted  winds,  especially  in  orbit  6681. 
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ORBIT  -  656S.1 

date  -  12286 


SZA(d«9)  51 .» 
GLAT(4«9)  -45.7 


41.7  44.6  58.6  81.2 

-19.6  11.5  38.2  68.9 


0«8«T  »  6681.0  SZA{<Je9)  46.0  46.6  52.0  60.8  71.8 

DATE  «  82294  GUT(<le9)  -24.5  -6.0  12.3  30.6  48.6 


OPeiT  =  6604.0  SZA(<leq)  44  0  45.4  47.3  57.9  73.2 

DATE  .  82239  GLAr(deg)  -29.3  -11.0  13  4  31  5  55.6 


figure  6.4  Three  orbits  of  DE2  showing  the  surface  plot  at  left  an  an  altitude  slice  at  110  km.  The 
altitude  slice  shows  the  VSH  model  meridional  winds  (plusses),  the  HWM  winds  (triangles)  and 
the  invert  green  line  winds  (boxes). 
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CHAPTER? 


ALTITUDINAL  PROFILES 


7.1  Red/Green  Line  Profile  Comparisons 

During  the  latter  orbits  of  the  DE2  mission  the  FPI  was  operated  in  a 
dual  mode  measuring  both  red  line  (6300  A)  emissions  and  green  line  (5577 
A)  emissions  during  each  second  of  operation.  The  red  line  winds  can  be 
taken  to  be  true  winds  at  each  altitude  with  no  inversion  procedure 
necessary.  In  this  chapter  these  red  line  higher  altitude  winds  are  stacked 
upon  lower  altitude  green  line  winds  to  test  for  consistency  and/or 
differences. 

Figure  7.1  shows  the  last  six  profiles  of  orbit  8072  with  the  red  line 
plotted  above  the  solid  line  at  250  km  and  the  inverted  green  line  winds 
below.  There  is  a  remarkable  correlation  in  the  wind  speeds  at  the  point 
where  they  meet  at  250  km.  Specifically,  profiles  28  and  29  almost  have 
identical  red  and  green  line  speeds  at  the  jimction  at  250  km  while  the  other 
orbits  match  the  wave  pattern  well.  This  orbit  is  during  a  low  period  of 
geomagnetic  activity  and  passes  through  high  latitudes  (-65°N  to  41°N). 
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figtire  7.1  Six  profiles  of  DE2  orbit  8072  showing  red  line  (6300A)  meridional  winds  above  250  km 
and  inverted  green  line  meridional  winds  below  250  km.  The  points  are  fitted  with  a  cubic  spline 
interpolation. 


figure  7.2  The  same  six  profiles  as  in  figure  7.1  plotted  sequentially.  The  two  strainght  lines 
show  the  downward  phase  propagation  in  the  wind  peaks  t3q>ical  of  gravity  waves. 
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Figure  7.2  stacks  these  profiles  next  to  each  other  sequentially.  The 
two  lines  drawn  through  the  peaks  in  the  profiles  show  a  typical  downward 
phase  propagation  characteristic  of  gravity  waves.  It  is  interesting  that 
this  is  evident  even  during  a  pass  through  auroral  latitudes  (since 
structure  tends  to  become  more  chaotic  there)  and  across  the  250  km 
red/green  line  boundary  in  these  plots. 

Figures  7.3  and  7.4  similarly  show  the  profiles  for  DE2  orbit  8092. 
This  orbit  passes  through  the  southern  auroral  region  (-85°  S  to  -62°  S).  The 
downward  phase  propagation  is  not  as  easy  to  identify  in  figure  7.4  for  this 
orbit  as  it  was  for  orbit  8072,  but  a  case  can  be  made  for  some  wave 
structure. 

An  attempt  was  made  to  locate  an  orbit  as  close  to  the  noon  to 
midnight  meridian  as  possible.  The  FPI  had  to  be  turned  off  when  too  close 
to  this  track  since  it  would  have  been  looking  directly  into  the  sim.  Orbit 
8088's  track  was  within  33°  of  this  meridian  and  its  profiles  are  shown  in 
figure  7.5.  Along  this  track  is  where  the  largest  wind  speeds  should  be 
detected,  but  here  the  largest  speeds  are  still  only  ~350  m/s  (profile  #23). 
Although  these  speeds  are  larger  than  other  orbits  presented  here,  speeds 
approaching  500  m/s  would  not  be  imreasonable.  Other  orbits  should  be 
studied  to  find  greater  speeds  and  to  study  the  associated  altitudinal  wave 
structures. 

Finally  Table  7-1  presents  a  summary  of  the  latitudinal  and 
altitudinal  wavelengths  found  in  the  40  plus  orbits  processed  in  this  study. 
These  wavelengths  agree  with  those  of  Nardi  (1991)  quite  well.  It  is 
interesting  to  note  that  for  both  the  altitudinal  and  latitudinal  waves,  the 
medium  Ap  orbits  showed  less  deviation  than  the  low  Ap  orbits.  The 
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reason  for  this  is  uncertain  and  may  be  only  a  function  of  the  peirticular 
orbits  selected  for  this  study. 


DE2-FPI  ORBIT  8002  RED/GREEN  LINE  WINDS 

Profile  #1  Profile  #2  Profile  #3 


figure  7.3  Six  profiles  of  OE2  orbit  8092  showing  red  line  (6300A)  meridionsi  winds  above  200  km 
and  inverted  green  line  winds  below  200  km.  The  points  are  fitted  with  a  cubic  spline 
interpolation. 
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figure  7.4  The  same  six  profiles  as  in  figure  7.3  plotted  sequentially.  The  two  straight  lines  show 
some  downward  propagation  of  wave  activity,  although  not  as  strongly  as  orbit  8072. 


DE2-FPI  ORBIT  8088  RED/6SEEN  LINE  WINDS 

Profile  #18  Profile  #19  Profile  #20 
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TABLE  7-1 


GREENLINE  WAVES 

(INVERTED  MERIDIONAL  WINDS) 
LATITUDINAL  WAVES: 

■  Wavelength  (degrees  Lat^  SandanlDeviatiop  (Degrees 


U>WAP(<2()): 

12.18 

3.88 

MEDAP(2(M0): 

10.61 

2.80 

HIGHAP(>40): 

13.51 

7.35 

OVERALLAVG: 

11.92 

4.48 

ALTITUDINAL  WAVES: 


LGWAP: 

72  km 

10km 

62-82  km 

MEDAR 

73  km 

6km 

67-79  km 

HIGHAP: 

72  km 

15  km 

57-97  km 

OVERALLAVG: 

72.3  km 

10.3  km 

62.3 -84.6  km 
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CHAPTERS 


DISCUSSION 

A  wind  inversion  technique  has  been  developed  to  recover  daytime 
lower  thermospheric  meridional  winds  from  satellite  measurements  of 
[01]  5577  A  emissions.  A  data  base  of  forty-four  DE-2  orbits  have  been 
processed  and  analyzed. 

An  onion  peel  brightness  inversion  was  used  to  obtain  an  inverted 
volume  emission  rate  from  the  limb-scanned  brightness  measurements 
derived  from  raw  spectrogram  data  of  the  DE-2  Fabry-Perot  Interferometer. 
These  inverted  voliime  emission  rates  were  then  used  to  apply  a  least 
squares  fit  to  chapmem  layers  to  obtain  an  inverted  meridional  wind  profile. 

These  inverted  winds  were  compared  to  the  Vector  Spherical 
Harmonic  (VSH)  Model  and  the  Horizontal  Wind  Model  (HWM). 
Comparisons  to  the  VSH  Model  revealed  a  good  correlation  at  low  and  mid 
latitudes  with  much  poorer  and  variable  cot.  elation  at  high  latitudes.  Even 
at  low  and  mid  latitudes  the  VSH  model  is  unable  to  duplicate  the  small 
scale  gravity  wave  effects  on  the  inverted  winds,  even  though  the  model 
does  incorporate  semi-diurnal  tides. 

The  comparisons  to  the  empirical  HWM  model  do  not  differ 
significantly  with  the  VSH  model  comparisons.  Once  again,  the  best 
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correlation  is  at  low  and  mid  latitudes  with  the  poorest  correlation  at  high 
latitudes. 

As  with  the  latitude  comparisons,  the  models  fit  the  data  best  at  low 
and  mid  geomagnetic  activity  indices.  For  high  geomagnetic  indices,  the 
correlation  is  poorer.  Interestingly  enough,  there  appears  to  be  a  slightly 
better  fit  at  medium  Ap/Kp  values  than  at  low  Ap/Kp  values.  The  reason 
for  this  is  xmcertain  and  should  be  investigated  further. 

The  altitudinal  structure  of  the  inverted  winds  was  also  examined. 
The  inverted  meridional  winds  revealed  much  more  wave  activity  than  did 
the  line-of-sight  winds.  Red  line  derived  winds  (6300  A  emissions)  were 
stacked  on  top  of  the  inverted  green  line  emission  (5577  A)  and  compared. 
A  remarkably  good  fit  occured  between  the  two  profiles  at  the  altitude  where 
they  met  (200-250  km).  This  is  an  encouraging  result  confirming  the 
strength  of  the  DE-2  FPI  data.  More  wave  activity  than  was  expected  was 
foimd  in  the  red  line  winds  above  250  km. 

When  stacked  sequentially,  the  profiles  show  the  stsmdard  downward 
propagation  of  gravity  waves.  This  result  appeared  even  across  the  meeting 
point  of  the  red  and  green  line  derived  winds. 

The  highest  wind  speeds  revealed  in  this  study  were  approximately 
350  m/s.  These  speeds  occured  when  the  satellite's  track  was  close  to  the 
noon  to  midnight  meridian.  Wind  speeds  as  high  as  500  m/s  may  appear, 
however,  in  further  studies  of  lower  thermospheric  winds  under  active 
geomagnetic  conditions. 

Finally,  a  statistical  calculation  of  the  latitudinal  and  altitudinal 
wavelengths  was  conducted.  The  results  confirmed  Nardi's  (1991)  initial 
resiilts  with  the  latitudinal  waves  ranging  from  10  to  14  degrees  of  latitude 
and  the  altitudinal  waves  averaging  70  to  75  km. 
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APPENDIX  A 


A  SEARCH  FOR  EVIDENCE  OF  IMPACT  IONIZATION 
IN  THE  ATMOSPHERIC  EXPLORER  UNIFIED  ABSTRACT 

DATABASE 


A.l  Introduction 

This  study  investigated  the  Atmospheric  Explorer  Satellite  (Unified 
Abstract)  database  for  evidence  of  impact  ionization  of  neutral  species  upon 
contact  with  the  satellite  surface.  This  search  is  based  on  a  paper  by 
Hanson,  Sanatani,  and  Hoffman  (1981)  that  addressed  ionization 
measurements  from  the  retarding  potential  analyzer  (RPA)  and  the 
magnetic  ion  mass  spectrometer  (MIMS)  on  the  Atmospheric  Explorer 
Satellites.  This  initial  investigation  detected  impact  ionization  or  “ion 
sputtering”  from  the  raw  (unaveraged  and  improcessed)  MIMS  and  RPA 
data.  Impact  ionization  occurs  when  a  neutral  atom,  molecule,  or  ion 
strikes  a  (satellite)  surface  with  sufficient  kinetic  energy  to  ionize  that 
particle  or  an  ambient  surface  particle.  Until  now,  it  has  been 
undetermined  whether  this  ionization  effect  made  the  transition  from  the 
raw  data  into  the  Atmospheric  Explorer  unified  abstract  database. 
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A.2  Earlier  Cionclusioiis 

Hanson,  Sanatani,  and  Hoffman  (1981)  reported  impact  ionization 
resulting  in  alkali  ions  Na+  and  K+  and  molecular  ions  NO+  and  O2+. 
These  ionizations  took  place  during  the  night  orbits  at  low  altitudes 
(approximately  150  km). 

The  Magnetic  Ion  Mass  Spectrometer 

The  Magnetic  Ion  Mass  Spectrometer  (MIMS)  measured 
concentrations  of  ions  from  1  to  64  atomic  mass  units  (amu).  The  entrance 
aperture  looked  radially  out  from  the  satellite’s  equator.  Three  slits  behind 
the  entrance  grid  separated  the  ions  according  to  low,  mid,  or  high  mass  in 
the  mass  ratijs  1:4:16.  These  three  mass  ranges  were  measured  by  the 
instrument  at  the  same  time  (1  to  4,  4  to  16,  and  16  to  64)  over  a  1  second 
interval.  An  alternate  mode  used  a  9  second  scan  period  to  and  extended 
the  mass  range  to  90  amu.  Figure  A-1  is  a  sketch  of  the  entrance  geometry 
of  the  MIMS  analyzer. 

A.2.2  Alkali  Ions 

In  their  paper,  Hanson  et  al.  [1981]  concluded  that  the  source  of  Na 
and  K  atoms  could  be  the  ocean  salt  spray  from  the  launch  sites  in  Florida 
and  California.  Several  factors  support  this  thesis.  First,  they  suspect  that 
the  MIMS  entrance  grid  itself  was  the  main  source  of  Na+  and  K+  since 
these  concentration  measurements  are  heavily  dependent  on  ram  angle. 
Ram  angle  is  the  angle  between  the  normal  to  the  entrance  grid  and  the 


57 


trajectory  of  the  incoming  particle.  Therefore  a  particle  entering  the 
instniment  at  a  90  degree  angle  to  the  entrance  grid  forms  a  zero  degree 
ram  angle  and  a  particle  entering  from  the  side  of  the  MIMS  forms  a  large 
ram  angle  near  90  degrees.  Na  and  K  could  have  infiltrated  into  the 
satellite  payload  upon  takeoff  and  these  deposits  would  have  been  upon  the 
MIMS  entrance  grid  wires.  Na  and  K  on  the  front  side  of  the  grid  wires 
may  have  been  ionized  from  collisions  of  neutral  O2  and  N2  and  sputtered 
away  from  the  MIMS  analyzer.  This  hypothesis  agrees  with  the  low 
measurements  of  Na+  and  K+  at  low  ram  angles. 

For  larger  ram  angles,  larger  measurements  of  Na+  and  K+  were 
detected.  This  may  be  due  to  Na  and  K  on  the  back  sides  of  the  wires 
becoming  ionized  and  sputtering  into  the  MIMS  analyzer.  The  ionized  Na 
and  K  from  the  front  of  the  wires  (low  ram  angles)  have  a  better  chance  of 
escape  from  the  MIMS  since  they  face  away  from  the  instrument  than 
those  on  the  backs  of  the  wires  facing  into  the  instrument. 

If  the  Na+  and  K+  are  from  ocean  salt  spray  deposits,  their  lifetimes 
were  not  infinite  once  the  satellite  was  in  orbit  and  removed  from  the 
source.  The  measurements  of  these  ions  over  the  lifetime  of  the  satellite 
confirms  this.  The  exponential  decay  time  constant  is  about  500  orbits 
according  to  the  law: 


I(n)Ml(0)exp[-n/500] 


ION  ^CAN  AM^UTUOeiark.  unlltl 
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figure  A-1  The  Magnetic  Ion  Mass  Spectrometer.  This  drawing  shows  the  entrartce 
grids  and  slits,  the  magnet,  collector  slits,  and  the  electron  multipliers. 
The  ion  trajectories  are  shown  as  dotted  lines 
(NASA,  1987) 


figure  A-2  The  peak  alkali 
ion  current  near 
1 50  km  as  a 
function  of  orbit 
number  (or  AE-C. 
(Hanson  et  al.,  1981) 
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where  1(0),  I(n)  are  the  initial  and  final  concentrations  lespectively 
and  n  is  the  orbit  number.  Figure  A-2  shows  a  plot  of  this  decay  for  Na+ 
and  K+. 

What  about  ambient  sources  of  Na  and  K  (sources  at  150  km  altitude)? 
These  sources  were  ruled  out  due  to  the  lack  of  ram  energy  in  the  mass 
spectrum  measurements  and  due  to  the  decay  of  Na  and  K  measurements 
with  time  and  orbit  number. 


A.2.3  Molecular  Ions 

The  molecular  ions  NO+  and  O2+  were  detected  on  the  mass  ion 
spectrum  from  the  analog  AE  data  with  a  shift  of  one  atomic  mass  unit 
away  from  the  ambient  NO+  and  O2+  due  to  their  lack  of  ram  energy.  O2+ 
displayed  a  decrease  by  a  factor  of  e  every  350  orbits  similar  to  the  alkali 
ions.  NO+  showed  a  less  time  dependent  decrease,  but  did  decrease  with 
time  also  (See  figure  A-3). 

The  main  source  of  O2+  is  from  collisions  of  O2  with  the  satellite 
surface.  The  main  source  of  NO+  is  probably  from  0+  formation  from 
collisions  of  atomic  oxygen  with  the  satellite.  Once  the  0+  is  formed  it 
combines  with  a  nitrogen  atom  from  the  satellite’s  surface.  Where  does 
this  surface  N  atom  come  from?  Possibly  from  the  use  of  hydrazine  as  an 
on  board  propellant. 

There  was  also  no  dependence  on  ram  angle  for  the  formation  of 
NO+  and  02+-  This  is  a  reasonable  result  since  the  sources,  ambient  0  and 
O2,  are  not  deposits  on  the  entrance  grid  wires  as  was  the  case  with  Na  and 
K. 
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A.2.4  Summary 

The  efficiency  of  the  impact  process  is  about  one  ion  per  10®  N2 
collisions  and  10®  O 2  collisions,  and  one  ion  per  3x10®  0+  collisions.  The 
formation  of  NO+  ions  is  from  collisions  of  O  with  the  spacecraft  surface  . 
The  abundance  of  NO+  is  highly  variable  due  to  the  extreme  variability  in 
the  supply  of  surface  nitrogen. 


A.3  EfEects  of  Impact  Ionization 

None  of  this  would  matter  if  the  impact  ionization  had  no  effect  on  the 
spacecraft  or  its  instruments.  The  currents  produced  by  these  ions, 
however,  did  cause  problems  in  the  operation  of  some  sensors  and  in  the 
interpretation  of  the  raw  data.  If  we  can  better  understand  and  anticipate 
these  currents,  we  can  design  instruments  to  correct  for  these  currents  or 
even  use  them  to  assist  in  the  data  gathering  process. 

A.4  The  Unified  Abstract  Atmospheric  Explorer  Data 

I  examined  the  MIMS  O2+  number  density  data  to  determine 
whether  any  of  the  sputtered  ions  were  evident  in  the  data.  Evidence  might 
possibly  have  shown  up  in  the  form  of  increased  number  density  at 
approximately  150  km  or  in  any  other  unusual  deviations  in  the  altitude 
profile  at  150  km. 

A  first  guess  hinted  that  all  impact  ionization  would  have  been 
filtered  out  of  the  UA  database  since  the  impact  ions  had  an  additional  ram 
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energy  that  separated  them  from  ambient  ions  by  one  atomic  mass  unit. 
This  indeed  proved  to  be  the  case.  The  raw  data,  however,  was  not  saved 
and  cannot  be  studied  further. 

The  MIMS  O2+  was  studied  from  the  AE-C  satellite.  The  satellite 
was  in  orbit  from  December  1973  to  December  1978.  Since  the  impact 
ionization  was  most  evident  in  the  first  part  of  the  satellite’s  life,  I 
examined  only  the  first  thousand  orbits. 

Figure  A-4  shows  a  crude  plot  of  O2+  number  density  vs.  orbit 
number.  No  data  was  available  until  Julian  day  74026  on  the  AE-C  MIMS 
O2+.  The  data  displayed  represents  altitudes  from  140  to  160  km.  The  gap 
around  orbit  711  to  875  is  another  period  where  no  data  was  available,  but 
the  average  values  should  be  decreasing. 

The  overall  decrease  is  not  a  reflection  of  the  impact  ionization 
decrease  as  discussed  above,  but  is  simply  a  reflection  of  the  night  side 
ionospheric  decrease  as  the  satellite  moves  from  da3d^ime  to  nighttime 
measurements.  Even  if  any  impact  ionized  ions  were  in  this  data, 
separating  their  decrease  with  orbit  number  from  the  decrease  due  to  the 
satellite’s  local  solar  time  would  be  difficult.  In  any  future  data,  sputtered 
ion  measurements  should  be  cataloged  and  stored  separately  from  ambient 
ions. 

A.5  Conclusions 

No  detectable  impact  ionization  evidence  exists  in  the  Unified 
Abstract  Atmospheric  Explorer  database.  Although  a  negative  conclusion, 
this  study  does  validate  the  reliability  of  the  MIMS  ion  data  as  a  true 
measure  of  the  E  region  ambient  ions.  This  could  be  an  important 


63 


consideration  in  future  models  of  the  lower  thermosphere  and  mesosphere 
that  would  include  ion  number  densities. 
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